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Abstract. Two-dimensional quadratic algebras are generalizations of Lie algebras that 
include the symmetry algebras of 2nd order superintegrable systems in 2 dimensions as spe¬ 
cial cases. The superintegrable systems are exactly solvable physical systems in classical 
and quantum mechanics. Distinct superintegrable systems and their quadratic algebras can 
be related by geometric contractions, induced by Inonu-Wigner type Lie algebra contrac¬ 
tions. These geometric contractions have important physical and geometric meanings, such 
as obtaining classical phenomena as limits of quantum phenomena as —>■ 0 and nonre- 
lativistic phenomena from special relativistic as c —>■ oo, and the derivation of the Askey 
scheme for obtaining all hypergeometric orthogonal polynomials as limits of Racah/Wilson 
polynomials. In this paper we show how to simplify the structure relations for abstract 
nondegenerate and degenerate quadratic algebras and their contractions. In earlier papers 
we have classified contractions of 2nd order superintegrable systems on constant curvature 
spaces and have shown that all results are derivable from free quadratic algebras contained 
in the enveloping algebras of the Lie algebras e(2, C) in flat space and o(3, C) on nonzero 
constant curvature spaces. The quadratic algebra contractions are induced by generaliza¬ 
tions of Inonii-Wigner contractions of these Lie algebras. As a special case we obtained the 
Askey scheme for hypergeometric orthogonal polynomials. After constant curvature spaces, 
the 4 Darboux spaces are the 2D manifolds admitting the most 2nd order Killing tensors. 
Here we complete this theoretical development for 2D superintegrable systems by showing 
that the Darboux superintegrable systems are also characterized by free quadratic algebras 
contained in the symmetry algebras of these spaces and that their contractions are also 
induced by Inonii-Wigner contractions. We present tables of the contraction results. 

Key words: contractions; quadratic algebras; superintegrable systems; Darboux spaces; 
Askey scheme 

2010 Mathematics Subject Classification: 22E70; 16G99; 37J35; 37K10; 33C45; 17B60 


1 Introduction 

We define an abstract nondegenerate (quantum) quadratic algebra as a noncommutative asso¬ 
ciative algebra generated by linearly independent operators H, Li, L 2 such that H is in the 

*This paper is a contribution to the Special Issue on Exact Solvability and Symmetry Avatars in honour of 
Luc Vinet. The full collection is available at http://www.emis.de/journals/SIGMA/ESSA2014.html 
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center, R = [Li, L 2 ] / 0 and the following relations hold: 

[Lj,R]= Y. e,>0, Ll = I, 

0<ei+e2+e3<2 

for some MA,e 2 ,e 3 £ C, and where [ 2 I, i?] = AB — BA is the commutator and {^ 1 ,^ 2 } = 
L 1 L 2 +L 2 L 1 is the symmetrizer. Also the operator R^ is contained in the algebra of symmetrized 
products: 

0 <ei+e 2 +e 3<3 

for some A^ei ,62,63 £ C. 

An abstract degenerate (quantum) quadratic algebra is a noncommutative associative algebra 
generated by linearly independent operators X, H, Li, L 2 such that H is in the center and the 
following relations hold: 

[X, L,] = p(Y,e3,e,LTL^2H^^X^\ J = 1, 2, 

0 < 6 i+e 2 +e 3 + 64 <l 

for some Rii,e 2 , 63,64 £ C. The commutator [Li,L 2 ] is expressed as 
[Li,L 2] = Qe4,62,63,64{^r^?,^}^'^^'"" 

0<6i+62+63+64<1 

for some <561,62,63,64 £ C. Finally, there is the relation: 

G= Y Se,,e2,e3,eALT,Ll\X^YH^"=^^ X^ = = I, 

0<61+62+63+64<2 

for some 5'ei,e2,63,64 £ C and where {L ^, is the 6-term symmetrizer of three operators. 

For both quantum quadratic algebras there is a natural grading such that the operators H, Lj 
are 2nd order, X is 1st order and 

ord([A, B]) < ord(A) -|- ord(iJ) — 1, ord(Ai?) = ord(A) -|- ord(R), (1.1) 

ord(/) = 0, ord(A + B) = max{ord(A), ord(iJ)}, ord(cA) = ord(A), 

for operators A, B, identity operator I and scalar c, with A ^ —B, A,B^0, and 0. Thus R 
is usually 3rd order, expression G is 4th order and F is 6th order. The field of scalars can be 
either M or C. 

There is an analogous quadratic algebra structure for Poisson algebras. An abstract nonde¬ 
generate (elassieal) quadratie algebra is a Poisson algebra with functionally independent genera¬ 
tors %, Cl, C 2 such that all generators are in involution with % and the following relations 
hold: 

{£„7^}= efc>0, £« = 1, 

0<6 i +62+63<2 

n^-F = n^- Y Xe,,e2,63^?^2"^"^ = 0 

0<6i+62+63<3 

for some Mij^e 2 , 63 , Xe^,e 2 ,e 3 £ C. An abstract degenerate (elassieal) quadratic algebra is a Poisson 
algebra with linearly independent generators X, %, C\, C 2 such that all generators are in 
involution with H and obey structure equations 

{x,c,} = Y ^i^62,63,64^r^r^'^^'^^^ 

0 < 6 l+ 62 -|-e 3 -|- 64 <l 


j = l,2 
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E 1^61,62, 63 ,64-^1 ^2 ^ 5 

0<ei+e2+e3+e4<l 

g= Y. = 0, = 1 (1.2) 

0< 61+62+63+64 <2 


( 7 ) 

for some Pei, 62 , 63 , 64 , < 361 , 62 , 63 , 64 ) 5 'ei, 62 ,e 3,64 £ C- There is a grading for these quadratic algebras 
with properties analogous to (1.1), but with the Poisson bracket instead of the commutator. 

These quadratic algebra structures arise naturally in the study of classical and quantum su¬ 
perintegrable systems in two dimensions and are key to the exact solvability of these systems, 
e.g., [2, 3, 6, 7, 9, 10, 25, 28, 32, 33]. A quantum 2D superintegrable system is an integrable 
Hamiltonian system on an 2-dimensional Riemannian/pseudo-Riemannian manifold with poten¬ 
tial that admits 3 algebraically independent partial differential operators commuting with H, 
the maximum possible. 

H = A + V, [H, Lj] = 0, = H, j = 1, 2,3. 

(In 2 dimensions we can always find Cartesian-like coordinates xi, X 2 such that 

and we adopt these coordinates in the following.) A system is of order k if the maximum order 
of the symmetry operators Lj (other than H) is k] all such systems are known for A; = 1, 2 [14, 
15, 16, 17, 23]. Superintegrability captures the properties of quantum Hamiltonian systems 
that allow the Schrodinger eigenvalue problem P'k = P'k to be solved exactly, analytically and 
algebraically. A classical 2D superintegrable system is an integrable Hamiltonian system on an 2- 
dimensional Riemannian/pseudo-Riemannian manifold with potential that admit 3 functionally 
independent phase space functions P, Ti, £2 in involution with P, the maximum possible: 

^ = + = Cs = n, j = 1,2,3, 

expressed in local Cartesian-like coordinates xi, X 2 , Pi, P 2 - A system is of order k if the maximum 
order of the constants of the motion Lj, j ^ 3, as polynomials in pi, p 2 is k. Again all such 
systems are known for k = 1,2, and there is a 1-1 relationship between classical and quantum 
2nd order 2D superintegrable systems [14, 15, 16, 17]. 

The possible superintegrable systems divide into four classes: 

1. First order systems. These are the (zero-potential) Laplace-Beltrami eigenvalue equations 
on constant curvature spaces. The symmetry algebras close under commutation to form 
the Lie algebras e(2,M), e(l,l), o(3,M) or o(2,l). Such systems have been studied in 
detail, using group theory methods, e.g., [27, 31]. 

2. Free triplets. These are superintegrable systems with zero potential and all generators 
of 2nd order. The possible spaces for which these systems can occur were classified by 
Koenigs [23]. They are: constant curvature spaces (6 linearly independent 2nd order sym¬ 
metries, 3 1st order), the four Darboux spaces (4 linearly independent 2nd order sym¬ 
metries, 1 1st order), and eleven 4-parameter Koenigs spaces (3 linearly independent 2nd 
order symmetries, 0 1st order). In most cases the symmetry operators will not generate 
a quadratic algebra, i.e., the algebra will not close. If the system generates a nondegenerate 
quadratic algebra we call it a free quadratic triplet. 





4 


R. Heinonen, E.G. Kalnins, W. Miller Jr. and E. Subag 


3. Nondegenerate systems [14, 15, 16, 17]. Here all symmetries are of 2nd order and the space 
of potentials is 4-dimensional: 

H(x) = aiV(;i)(x) -Fa 2 V( 2 )(x) -Fa 3 E( 3 )(x) -^ 04 . 

The symmetry operators generate a nondegenerate quadratic algebra with parameters aj. 

4. Degenerate systems [18]. There are 4 generators: one 1st order X and 3 second order H, 
Li, L 2 . Here, X'^ is not contained in the span of H, Li, L 2 . The space of potentials is 
2-dimensional: H(x) = aiV(i)(x) -|- 02 . The symmetry operators generate a degenerate 
quadratic algebra with parameters aj. Relation (1.2) is an expression of the fact that 4 
symmetry operators cannot be algebraically independent. 

Every degenerate superintegrable system occurs as a restriction of the 3-parameter poten¬ 
tials (i.e., 4-dimensional potential space) to 1-parameter ones, such that one of the symmetries 
becomes a perfect square: L = X^. Here X is a first order symmetry and a new 2nd order 
symmetry appears so that this restriction admits more symmetries than the original system. 

Strictly speaking, a nondegenerate 2D superintegrable system, both classical and quantum, 
is not a single system but in fact a family of superintegrable systems parameterized by three 
parameters, oi, 02 , 03 . Similarly a degenerate 2D superintegrable system, both classical and 
quantum, is a family of superintegrable systems parameterized by one parameter, 04 . 

Eor a quadratic algebra that comes from a nondegenerate 2D superintegrable system (classical 
and quantum) the constants Mei]e 2 ,e 3 and Xej^e 2 ,e 3 are polynomials in the parameters ai, 02 , 03 
of degree 2 — 61 — 62 — 63 and 3 — 61 — 62 — 63 , respectively. If all parameters aj = 0 the 
algebra is free. For a quadratic algebra that comes from a degenerate 2D superintegrable system 
(classical and quantum) the constants Rii,e 2 ,e 3 ,e 4 , Qei,e 2 ,e 3 ,e 4 and 5 e 4 ,e 2 ,e 3 ,e 4 are polynomials in ai 
of degrees 1 — 61 — 62 — 63 — 64 , 1 — 61 — 62 — 63 — 64 and 2 — 61—62 — 63 — 64 , respectively. If 
all parameters aj = 0 these algebras are free. 

Basic results that relate these superintegrable systems are the closure theorems: 

Theorem 1.1. A free triplet, classical or quantum, extends to a superintegrable system with 
potential if and only if it generates a free quadratic algebra Q. 

Theorem 1.2. A superintegrable system, degenerate or nondegenerate, classical or quantum, 
with quadratic algebra Q, is uniquely determined by its free quadratic algebra Q. 

These theorems were proved in [21], except for systems on the Darboux spaces which will 
be proved in this paper. The proofs are constructive: Given a free quadratic algebra Q one 
can compute the potential V and the symmetries of the quadratic algebra Q. Thus as far as 
superintegrable systems on specific spaces are concerned, all information about the systems is 
contained in the free quadratic algebras. 

We will refer to quadratic algebras associated with superintegrable systems as geometric, 
a subset of abstract quadratic algebras. For quadratic algebras associated with quantum super¬ 
integrable systems the order of A is its order as a partial differential operator. For quadratic 
algebras associated with classical superintegrable systems the order of A is its order as a poly¬ 
nomial in the momenta. Although there is a 1-1 relationship between classical and quantum 
geometric systems the corresponding classical and quantum geometric quadratic algebras are 
not the same; only the highest order terms in the structure equations agree. 

2 Contractions 

The notion of contractions for quadratic algebras is based on that for Lie algebras, e.g., [13, 29, 
34]: 
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Definition 2.1. Let g be a complex Lie algebra with an underlying vector space V and Lie 
brackets [ , ]• In the following we simply write it as g = (R, [ , ]). Suppose that for any 
e G (0,1], te: R —)• R is a a linear invertible operator and that lim converges for 

E— 

any X,Y G V. We use the notation 

lim t-%X,t,Y] = [X,R]o. 

€^0 + 

Then [ , ]o are in fact Lie brackets on R and we denote this Lie algebra by go = (R, [ , ]o). We 
say that go is a contraction of q (that is realized by the family of linear maps {te}eg(o,i]) and 
we denote it by g —)• go- 

Thus, as e —0 the 1-parameter family of basis transformations can become singular but the 
structure constants go to a finite limit. 

Note. In this paper all of the Lie algebra contractions needed are of standard Inonii-Wigner 
type [29], but for our theory on contractions of quadratic algebras and second order superinte¬ 
grable systems on constant curvature spaces [21] IW-contractions are not enough. In most cases 
we employ a generalized Inonii-Wigner contraction (Doebner-Melsheimer type) [4, 29], but in 
some specific cases we are forced to use a general contraction of Lie algebras in the sense of 
Saletan [30] . In a paper under preparation we shall demonstrate that all of these generalizations 
of Inonii-Wigner contractions are induced by a well defined family of contractions of the confor¬ 
mal Lie algebra so(4, C) to itself that follow from limiting processes for R-separable coordinate 
systems for wave equations, introduced by Bocher in his famous 1894 thesis [I]. We call these 
Bocher contractions. 

Of particular interest to us are contractions that are induced by e-dependent local analytic 
coordinate transformations Xj{e, x[,..., x'^f), j = 1,..., n on a manifold A4 such that the Ja¬ 
cobian det(^) 7 ^ 0 for e G (0,1], but the Jacobian is undefined or nonsingular in the limit as 
e ^ 0. 

Definition 2.2. Suppose that g = (R, [ , ]o) is a contraction of g = (R, [ , ]) that is realized by 
the family of linear maps {te}eg(o,i] • Let AI be a smooth manifold with a local coordinate system 
{xi,X 2 , ■ ■ ■ ,Xn) and let ■0: g —>■ C°°{TM.) be an embedding of Lie algebras, where C°°{TM) is 
the space of smooth functions on the cotangent bundle (the phase space) of M. which is equipped 
with its canonical Poisson brackets as Lie brackets. Suppose that Xj{e, x'^, ..., x'^f), j = 1,..., n 
are e-dependent local analytic coordinate transformations such that the Jacobian det(^) / 0 
for e G (0,1]. If the limit 

lim il){t^{v)){xi{e,x'^, .. .,x'J,.. .,Xn{e,x\,.. .,x'J) 
e—>0+ 

converges for any v G V and defines a Lie algebra homomorphism from go into a local expression 
of a certain space of {'TXi') for some smooth manifold Ai'. Then we say that the contraction 
g —)• go is implemented by Xj{e, x'l, ..., x'„), j = 1,... ,n and call this procedure a geometric Lie 
algebra contraction. This distinction between abstract contractions of Lie algebras and geometric 
contractions has been recognized from the earliest days of the theory, e.g., [13]. 

We give some pertinent examples. 

Example 2.1. Consider the complex three dimensional Lie algebra G3 defined by basis elements 
{'Pi,T’ 2 , T’} that satisfy [Pi, P 2 ] = 0, [Pi, P] = Pi, [P 2 , P] = P 2 - This algebra admits an abstract 
Inonii-Wigner contraction defined by te(Pi) = ePi, te(P 2 ) = eP 2 , L(P) = P- In this case the 
contracted Lie algebra, ^3o, coincides with Q3 and [Pi,P 2 ]o = 0, [Pi,P]o = Pi, [P 2 ,P]o = P 2 - 
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Now considering the complex analytic manifold with coordinates (xi,X 2 ) we can reali¬ 
ze ^3 by 


'll){Vl){xi,X2,Pxi,Vx2) =Pxi, 'lp{V2){xi,X2,Pxi,Px2) = Px 2 , 
i}{V){xi,X2.,Pxi,Px2) = XlPxi +X2Px2- 


A geometric implementation of the contraction is obtained by the substitution xi = ex'^, X 2 = 
ex' 2 - Then, taking the limit, we find 

Y\ui i){t^{Vi)){ex'i,ex'2,p^x\,Pex')= 1™ ePcx' = 1™ =Px\^ 

e-5>0+ 1 ^ e-5>0+ 1 e-5>0+ 6 ^ l 

Y\uijl){t^{V 2 )){ex'i,ex' 2 ,p^x\,Pex') = = lim , 

e-5>0+ 1 2 e-5>0+ 2 e-5>0+ 6 ^ 2 

Ym. i;{t^{V)){ex\,ex'2,Pex\,Pex') = lini €a;'i + ^A^Px' = x'lPxL + APx'y 

Though this Lie algebra contraction acts like the identity map here, we shall see that its action 
on Darboux quadratic algebras is nontrivial. 


Example 2.2. We again consider the Lie algebra ^3 with the same geometric identification. 
This algebra admits another abstract Inonii-Wigner contraction defined by teiVi) = 'Pi, 
L(P 2 ) = P 2 , A(P) = eP- The contracted Lie algebra, Q^o, is given by [Pi, P 2 ]o = 0, [Pi, P]o = 0, 
[P 2 ,P]o = 0. A geometric implementation of the contraction is obtained by the substitution 
xi = x'i + \ , X 2 = x' 2 - Then, taking the limit, we find 


Y\m. 'il){te{Vi)){x'i + -,x' 2 ,p^: ,i,Px') = lim p , , 1 = lYm p^i =Px',, 
e- 5 > 0 + e 2 ^_^Q+ ^ 1 -r^ ^_^Q+ 1 1 

lim V’ (A (P 2 )) (a; 1 + -, a :^2 > Px\+i^Px') = lini Px' = Px', 

e- 5 > 0 + e 2 ^_^Q+ 2 2 

^hm V’(A(P))(P2) (x[ + X2,P3,/+1 ,Pa;^^ = e T (^x[ + X2Px'^^ = Px[ ■ 


Hence the resulting map from QSq that is given by Pi 1 —)■ p^'^, V 2 1 —> Px^; Pi Px[ is a Lie 
algebra homomorphism with kernel spanned by Pi — P. Though this geometric implementation 
is not an isomorphism of the contracted Lie algebra we shall see that its action on Darboux 
quadratic algebras is isomorphic. 


There are exactly analogous implementations of geometric contractions in the quantum case. 

Definition 2.3 (algebraic contraction of quadratic algebras). Let A be a classical nondegenerate 
quadratic algebra with a generating set which consist of a Hamiltonian, P (second order element 
which lies in the center of A) and two second order constants of motion £ 1 , £ 2 - Let {£ 1 , £ 2 }^ = 
= P(P,£i,£ 2 ) be the Casimir of A. For any e G (0,1] consider a matrix G GL(3,C) of 
the form 

/2lqi(e) Ai^2(e) ^i,3(e)\ 

Ae = j A 2 q(e) A 2 , 2 (e) ^ 2 , 3 (e) 1 £ GL(3). 

V 0 0 As,3(e)/ 

Assume that the map from (0,1] to GL(3,C) that is given by e 1 —)■ A^ is continuous. For any 
e G (0,1] we have another set of generators, {£f, £|,P'^}, for A that is defined by 
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and satisfies the Casimir relation 


{C\,C\Y = = nYAnie)A22(e) - AMA^iY)) 


-2 


— (^ii(e)^22(e) — ^i2(e)^2i(e)) A 




Y. a,M^mf(ciy{n‘Y 

i+j+k=3, 0<ij,k 


If lim aijk{e) exists for any i, j, k we denote it by aj,fc(0). Then there exists a quadratic 
£->0+ 

algebra, with a set of generators {C\, that satisfy 

[ClClY = ( 7 ^ 0 )' = 

i+j-|-/c=3, 0<i,j,k 


we call ^0 the contraction of A with respect to {^e}ee(o,i]- 

Note that we can expand as quadratic expressions in C\, £|, For 

a contraction it might seem that we must also require these expansion coefficients to have 
finite limits as e —)• 0. However from the results of Section 3.1 the convergence of these other 
structure equations follows from the convergence of the Casimir. There is a completely analogous 
definition of contraction for quantum quadratic algebras. 

Just as for abstract classical and quantum Lie algebra contractions there are abstract classical 
and quantum quadratic algebra contractions that are induced by e-dependent local analytic co¬ 
ordinate transformations Xj(e, x'l, ..., x^), j = 1,..., n on a manifold Ai such that the Jacobian 
det(^) A 0 for e G (0,1], but the Jacobian is undefined or nonsingular in the limit as e —>■ 0. 
If an algebraic contraction A ^ B can be implemented by some coordinate transformation, we 
say that it is a geometric quadratic algebra contraction. In fact, all of the quadratic algebra 
contractions for Darboux systems discussed in this paper are geometric implementations. We 
will give many examples in the following sections. 

The notion of contraction applied to structures other than Lie algebras is not new, see for 
example [11] (and references there-in) and [26]. 

In [21] Lie algebra and quadratic algebra contractions for superintegrable systems on constant 
curvature spaces were related: 

Theorem 2.1. Every Lie algebra contraction of Q = e(2,C) or Q = o(3,C) induces a geometric 
contraction of a free geometric quadratic algebra Q based on Q, which in turn induces uniquely 
a contraction of the quadratic algebra Q with potential. This is true for both classical and 
quantum algebras. 


Here we will demonstrate the analogous result for Darboux spaces, using the conformal 
symmetry algebra G3 with basis {dx, dy, xdx + ydy}. 


3 Structure relations 


Although the full sets of classical structure equations can be rather complicated, the function E 
contains all of the structure information for nondegenerate systems and Q (only unique up to 
a nonzero scalar multiple) most of the information for degenerate systems. In particular, it is 
easy to show that [8, 21], {Ci,TZ} = {C 2 ,'R.} = so the Casimir contains within 

itself all of the structure equations. For degenerate systems we have [21] 


{X,Ci} = K 


dQ 

dLY 


{X,L2} = -K 


dQ 

dCY 


{A,A} = A-^. 
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where {H^K} = 0 . Here, R is a scalar, unless {X,Ci} and {X,C2} are linearly dependent. In 
the latter case there would exist 3 algebraically independent elements of the algebra in involution, 
including the Hamiltonian. This is impossible for a Hamiltonian system. Thus, except for some 
abstract quadratic algebras unrelated to geometric superintegrable systems, K will always be 
a scalar that can be normalized to 1. 

3.1 The quantum operators F and G 

The quantum case is similar to the classical case, but more complicated. From the Casimir 
relation 

R^-F = R^- (biL\ + b 2 Ll + + 64 ( 1 ^?, L2} + 65 ( 1 ^ 1 , 1 ^ 2 } + hLiL 2 Li + 67 ^ 2 ^ 11^2 

+ bfiH{Li,L2} + b^HLl + bioHLl + buH^h + bi2H^L2 + b^Lf + buL^ 

+ bib{Li, L2} + bieHLi + bnHL2 + + &i 9 -bi + ^20-^2 + b2iH + 622) = 0, 

we want to determine the structure relations. Noting that R = [Li,L2] and using operator 
identities 

L2RL2 = \[[L2,R],L2] + \{R,Ll], LiRLi = l[[Li,R],Li] + \{R,LI], 

L1RL2 + L2RL1 = -1[L2, [Li, R]] - i[Li, [L2, R]] + \{R, {Li, L2}}, ( 3 . 1 ) 

[Li,R 2 ] ={R,[Li,R]}, [L 2 ,R 2 ] ={R,[L2,R]}, 

[LUL2L1L2] = {[Li,L2],HLi,L2}}, 

XL1L2 + L2L1X = 1 {X, {Li, L2}} + l[X, [Li, L2]], 

(true formally for all operators Li, L2, R, not just for R the commutator), and setting 

[Li, R] = AiL\ + A2L2 + A^H^ + H4{Li, L2} 

+ A^HLi + AqH L2 + AjLi + AqL2 + AgH + Aiq, 

[L2, R] = RiL? + B2LI + + R 4 {Li, L2} 

+ B^HLi + BqHL2 + BjLi + BsLg + BgH + Rio, 

we can write [Lj, R^ ~ F] = 0 in the form {R, Qj} = 0 for some explicit differential operator Qj. 
This can only hold if Qj = 0 . As a result we find: 

Hi = 64 + Ifte, A2 = §62, A3 = 5612, A4 = 65 + ^67, 

A5 = bg, Aq = 610, A^ = 615 — I6162 + 64^5 — \b&bj, 

As = 614 — 5 & 2 fe 4 + bl — bgbe + 56567, Ag = 5617 — 36269 + 56568 — 360610, 

Alo = 5620 — 462613 — 3616265 

- 46264 + 565615 + 5646I + 4616267 + 36260 - 36560 - 4656067 - 460614, 

Ri = —561, R2 = —65 — 567, R3 = —5611, R4 = —64 — 560, 

R5 = —69, Ro = — bg , R7 = —613 + 56165 — 6| + 6167 — 56460, 

Bg = —615 + 36162 — 6465 + 36067, Bg = —5610 + 361610 — 56463 + 36769, 

Rio = —5619 + 461614 + 3616264 + 3616I — 564615 

— 56465 — 3616260 — 36167 + 46467 + 3646067 + 367613. 

For quantum degenerate systems, in the Casimir relation 
G(Li,L2,R,A,a) = 0 
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for 2 nd order superintegrable systems with degenerate potentials we assume that G is given, up 
to a multiplicative factor, and set 

G = ciL\ + C2L2 “ 1 “ “ 1 “ C 4 {-Li, ^2} “ 1 “ c^HLi + cqHL2 + } 

+ Cgjx^, L2} + CiqHX^ + CiiXLiX + Ci2^-h2^ + C13L1 + C14L2 
+ C15H + CiqX^ + Ci7, 

[X, Li] = CiLi + G2L2 + Ggff + C4X^ + C5, 

[X, L2] = DiLi + D2L2 + Dgff + D 4 X^ + D5, 

[Li, L2] = ^i{Li, X} + ^2{-h2, X} + E3HX + E4X^ + E5X. 

Using identities ( 3 . 1 ) we find 

[X,G] = {[X,Li],^i - ^{Gf + G2D4) - + D1D2) + 

+ {[X,L2],^2 + ^(^^iC'2 + C'2D2)-^(C2Di + d 2 ) + ^} = 0 , ( 3 . 2 ) 

[Li, G] = {[Li, L2], ^2 + j{CiC2 + G2D2) - |(C2 + C1D2) + ^(GiD2 - C2D1) 

+ ^ } + {[Li, X], ^3 + fciX} + {[L 2 , X], fc2X} = 0, (3.3) 

[L2, G] = {[Li, L2], -^1 + I(G1D2 + Dl) - j(GiDi + D4D2) + ^(G2Di - G1D2) 

- ^ } + {[Li, X], kgX} + {[L2, X], 2I3 + A:4X} = 0, ( 3 . 4 ) 

where 

Ai = ciLi + C 4 L 2 + + ^^8 H—X^, 

A2 = C2L2 + C4L1 + + ^cg H — X^, 

^3 = 2 c 7 X 3 + (c 8 + {Li,X} + (eg + {L2,X} + cioHX, 

ki = -cj{Cl + C2D4) + y (2G1G4 + G2D4 - G2E4) + ^(2G4Di + Gi^li - G1D4) 

+ ^(GiEi + G1D4 + D4E2 - G4D4) + Ci6, 
k2 = -C7(GiG2 + C2D2) + y (G2G4 - G2E2) + y (2G4D2 + G1G4 + Gi^;2) 

+ ^(G2D4 + G2E4 + D2E2 - G4D2), 

kg = -C7(GiDi + D1D2) + y (2G1D4 + D2D4 - D4E4) + ^{DiEi + D4D4) 

+ ^(G4Di - G1D4 - G1.B1 - D1.B2), 

k 4 = -C7(G2Di + Dl) + y (2G2D4 - G4D2 - D2E2) + y (2D2D4 + G4D1 + D4E2) 

+ ^{G4D2 - G2D4 - G2E4 - D2E2) + Ci6. 

Equating the coefficients of the 4 th order terms in ( 3 . 2 ) and the coefficients of the 5 th order 
terms in ( 3 . 3 ) and ( 3 . 4 ) we find 

[X,Li] =X^2 + G5, [X,L2] = -KA4+Dg, [14,12] = KAg + E^X, 

G 4 =Kc 4 , G2 = Kc2, C'3 = x|, C 4 = x(c9 + ^), 



10 


R. Heinonen, E.G. Kalnins, W. Miller Jr. and E. Subag 


Di = -Kcu D2 = -Kci, D 3 = -K^, E »4 =(c8 + ^) , 

Ei = J^(c8 + ^), E 2 = i^(c9 + ^), E3 = Kcio, E^ = 2 Kc7, 

for some constant K. Now we substitute these values back into (3.2), (3.3), and (3.4). We 
immediately see that k 2 = = 0 and 

ki = kA = ci 6 — 07(04 — 01 C 2 ) + Rr^(o 409 — 02 O 8 ) ^os + 

+ ^^(0408 — 0109) ^09 + 

and we obtain 

C 5 A 1 + D^A 2 + (^ - - 0102 )) {KA 2 + C 5 ) 

+ (t “ 

2 - J^2^(ol - 0102 )) {KA3 + E5X) - 2C5{A3 + kiX) 

+ {E5-Kki){X,A2} = 0, (3.6) 

-2 - K^^-f{cl - 0102 )) {KAs + E 5 X) - 2 D 3 {A 3 + k^X) 

-{E 5 -Kki){X,A^} = 0. (3.7) 

These equations split into terms of order 3,2,1 and 0. From equation (3.5) we find 

C5 = ^K-^K^{cl-ciC 2 ), D5 = -^K+^K^{cl-ciC 2 ) (3.8) 

except, possibly, for some degenerate cases. The condition that (3.8) is the unique solution 
of (3.5) is exactly that the set ([X, Ri], [X, S' 2 ]) is linearly independent. Otherwise the solution, 
though it always exists, is not unique. 

Substituting (3.8) into (3.6) and (3.7), we have 

2 C 5 {E 5 ^Kki) ^ ^ ^ 

2D,{E,^Kk,) ^ _ ^ 

whence we find 

E5 = Kciq — K^C7[c^ — C1C2) + Rr^(c4C9 — 0203) ^08 + 

+ K ^{ c 4 Cs - C1C9) (^Cg + . 

We conclude in both the classical and quantum cases that the Casimirs of superintegrable 
systems determine the structure equations. 

4 Free 2D 2nd order superintegrable systems 

As was shown in [14, 15, 16, 17] the “free” 2nd order superintegrable system obtained by setting 
all the parameters in a nondegenerate potential equal to zero retains all of the information 
needed to reconstruct the potential. Thus we can, in principle, restrict our attention to free 
systems. First we review from [14, 15, 16, 17, 21] how the structure equations for 2D 2nd 
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order nondegenerate classical superintegrable systems are determined. Such a system admits 

a symmetry £ = '^a^^PiPj + W if and only if {'H,C} = 0, i.e., the Killing equations are 
2 

satisfied and Vhi = A ^ Here Wi = dx^W with a similar convention for subscripts on V. 

i=i 

The equations for W can be solved provided the Bertrand-Darboux equation dx^W 2 = dx 2 Wi 
holds. For a superintegrable system with independent symmetries £i = 

C 2 = Eb^^PiPj + W^^\ we can solve the two independent Bertrand-Darboux equations for the 
potential to obtain the canonical system 


K 22 - Fii = Fi2 = A^'^Vi + B^'^V2 


,22 


12 i 


? 12 i 


(4.1) 


Here, 


= _G 2 + 


D 


( 2 ) 


D 


= -Gi - 


D 


( 0 ) 


D 


D = det 


6 ^^ — h 


^22 = 2 Gi + 


D 


( 3 ) 


D 


5^2 = - 2 G 2 - 


D 


( 1 ) 


all _ ^22 
,22 


ai 2 

612 


= det 


D( 3 ) = det 


12 


all _ „22 


3ai2 

a 2 


all _ ^22 


36}2 511 - 6 


D 


D(o) = det 


12 


3a. 
3612 


3a 

3612 511 _ 522 D 




-ai 2 
-612 
12 


D( 2 ) = det 


1 


12 


V3612 612 ;’ 


l 22 ) 


where A = exp G. If the integrability equations for (4.1) are satisfied identically then the solution 

3 

space is 4-dimensional and we can always express the solution in the form V (x) = (^) + 

i=i 

04 , where 04 is a trivial additive constant. In this case the potential is nondegenerate and 3- 
parameter. Another possibility is that the solution space is 2-dimensional with general solution 
K(x) = aiV(i)(x) -|- 02 - For nondegenerate superintegrability, the integrability conditions for 
the canonical equations must be satisfied identically, so that V, Vi, V 2 , Vu can be prescribed 
arbitrarily at a fixed regular point. 

To obtain the integrability conditions for equations (4.1) we introduce the dependent variables 
B/(i) = Fi, VF( 2 ) ^ ^ matrices 


/ 1 T(i)\ / 0 0 1 \ 

w = 1T(2) ^ ^( 1 ) ^ ^12 ^12 Q 

\it( 3 )/ \A 13 b 13 B ^^- A^y 

yll 3 = - Af + ^12^22 + ^12^12 - ^22^12^ 

B 13 = 5I2 _ ^22 ^ ^12^12^ ^23 ^ ^12 ^ ^12^12^ 


/Ai 2 B 12 0 \ 

A( 2 ) = A 22 B 22 1 1 

\A23 B23 A^y 


B23 ^ ^12 ^12^12 


Then the integrability conditions for system j = 1, 2, must hold: 


a(2) - 41) = a(i)a( 2) - A(2)yl{i) = 4(1), a(2)]. 


(4.2) 


If and only if (4.2) holds, the system has a 4D vector space of solutions V. 

There is a similar analysis for a “free” 2nd order superintegrable system obtained by setting 
the parameter in a degenerate potential equal to zero, [18]: The free system retains all of the 
information needed to reconstruct the potential. All such degenerate superintegrable systems 
with potential are restrictions of nondegenerate systems obtained by restricting the parameters 
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so that one 2nd order symmetry becomes a perfect square, e.g., L\ = . Then T" is a 1st order 

constant, necessarily of the form X = ^ipi + ^ 2 P 2 , without a function term. Since the degenerate 
systems are obtained by restriction, the potential function must satisfy the equations (4.1) 
inherited from the nondegenerate system, with the same functions In addition the 

relation {X^T-L} = 0 imposes the condition + ^ 2 X 2 = 0. By relabeling the coordinates, 

we can always assume ^2 7 ^ 0 and write the system of equations for the potential in the form 
E 2 = CVi, V 22 = Ell + V 12 = where 

C^(xi,X2) =C^^(xi,X2) = C^^(xi,X2) = 

42 42 42 

To find integrability conditions for these equations we introduce matrices 


V = 



= 


0 1 
0 92^2 + _ (^22 


= 


0 

0 


Then integrability conditions for system dxjV = B^-^^v, j = 1, 2, must hold: 


Rf ^ - B^^^ = ^(2)] _ 


(4.3) 


If and only if (4.3) holds, the system has a 2D space of solutions V. Since V = constant is 
always a solution, (4.3) is necessary and sufficient for the existence of a nonzero 1-parameter 
potential system. In this case we can prescribe the values V, V 2 at any regular point xq; there 
will exist a unique E(x) taking these values. 


4.1 Free triplets 

Here we review information about free triplets that was presented and proved in [ 21 ]. A 2nd 

order classical free triplet is a 2D system without potential, T-Lq = and with a basis 

2 .. 

of 3 functionally independent second-order constants of the motion ^ CL'^g'^PiPj-, = 

*j=i 

o|\, s = 1,2,3, £( 3 ) = Bq. Since the duals of these constants of the motion are 2nd order 
Killing tensors, the spaces associated with free triplets can be characterized as 2D manifolds 
that admit 3 functionally independent 2nd order Killing tensors. As mentioned above, they have 
been classified in [23]. Since the vectors {h(s)}, h(s)*’’(x, y, z) = (O(j^)) ^(s)) form a linearly 

independent set, there exist unique 3x3 matrices guch that j, s = 1 , 2 . 

2 

By linearity, any element £ = X] <A^PiPj of fh® space of 2nd order symmetries spanned by the 
basis triplet is characterized by matrix equations 

dxX = C^A\y^ j = l,2, h*’'(a:,y, 2 ;) = (a^^a^^a^^) . (4.4) 

In particular, at any regular point xq we can arbitrarily choose the value of the 3-vector hg 
and solve (4.4) to find the unique symmetry £ of Bq such that h(xo) = hg. A normalization 
condition for the (4.4) is valid for = 1/A, = 0, i.e., for Bq. From this and the 

requirement that the £ are free constants of the motion we find 



( -Gi 

-G 2 

0 \ 


(d? 

^(2) 

012 

-G2-£j?\ 

cW = 

iMA 

2'-'n 

1 1^(2) 

2CT1 2012 

1^(2) 

2*-'ll 

, £(2) ^ 

^(2) 

021 

^(2) 

O22 

^(2) 

O2I 


- 2Cg) 

9^(2) 
Ct2 ^(-22 

2C2?j 


'v 0 

-Gi 

-G2 / 
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with the 4 functions ^ > ^ 22 ^ define the functions A^'^, by 

the requirement 

d? = -1 g, - is 

then equations (4.4) agree with the equations that are obtained from a superintegrable system 
with nondegenerate potential satisfying (4.1). Thus, for a free system there always exist unique 
functions A^^, B^B Then necessary and sufficient conditions for extension to a system with 
nondegenerate potential V satisfying equations (4.1) are that conditions (4.2) hold identically. 

This analysis also extends, via restriction, to superintegrable systems with degenerate po¬ 
tential. A free triplet that corresponds to a degenerate superintegrable system is one that 
corresponds to a nondegenerate system but such that one of the free generators can be chosen 
as a perfect square. For these systems conditions (4.3) for the potential are satisfied identically. 

Similarly, we define a 2nd order quantum free triplet as a 2D quantum system without po¬ 
tential, Hq = y^(9ii -|- 522 ), and with a basis of 3 algebraically independent second-order 
symmetry operators 

1 2 

= diiXai^k^dj){yL), k = 1,2, 3, L 3 = i^o- 

* 4=1 

There is a 1-1 relationship between classical and quantum free triplets. 


Cij = ^Gi - ^A22, 
^2 3 3 

eg) = --G2- -B 
22 3 3 


5 Free Darboux systems 

The Darboux spaces admit a 1-dimensional space of Killing vectors and a 4-dimensional space of 
2nd order Killing tensors. Thus each space can admit at most one superintegrable system with 
degenerate potential, and each space does so. We merely need to check that equations (4.2) are 
satisfied. Then we can compute the degenerate potential. Turning off the 1-parameter potential 
produces a single free degenerate quadratic algebra which we list below. There are no more pos¬ 
sibilities. There are a number of possibilities for free triplets to define a nondegenerate quadratic 
algebra for a Darboux space, however. We classify the possibilities up to conjugacy under the 
action of the 1-parameter symmetry group of the manifold. Note that the 4-dimensional space 
of free constants of the motion is not obtained from the enveloping algebra of an underlying 
symmetry group. We shall see that there is a 1-1 relationship between free quadratic algebras 
and restrictions of quadratic algebras of nondegenerate superintegrable systems. We adopt the 
labeling of superintegrable systems and constants of the motion on Darboux spaces introduced 
in [19], using a tilde to differentiate between a free triplet and its associated superintegrable 
system. In the following sections, with the aid of MAPLE, we classify all possible free quadratic 
algebras generated by the 2nd order Killing tensors, up to conjugacy. Then, using MAPLE 
again, we verify for each quadratic algebra that the integrability conditions (4.3) are satisfied 
and we compute the nondegenerate potentials. Most of the results are presented in lists but in 
Section 5.2 we give more details on the construction of the superintegrable system with potential 
in one case. 

Each of the Darboux spaces can be embedded as a surface in 3 dimensions if we regard the 
ignorable variable as an angle, i.e., X = f{x) cos(y), Y = f{x) sin(?/), Z = h{x), and this is not 
unique [5]. We give an illustrative example for each case. 
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5.0.1 Free Darboux 1 systems 

The space Darboux 1 (Di) has free degenerate Hamiltonian 

DID: n=^{pl+pl) 

with a single Killing vector 1C = py and a basis, {T-L,K?, Xi, X 2 } for the 4-dimensional space of 
2nd order Killing tensors. Here, 

2 

Xl = PxPy - pI) > = PyivPx - xpy) -^{pI+ pI) ■ 

The commutation relations are 


{lC,Xi] = -2rL, {lC,X2] = Xi, {T’pT’s} =-2/C^ 

and there is the functional relation XKX 2 + X^ -|- = 0. The degenerate potential is V(x, y) = 

^ + ^ 2 - 

As shown in [20] , a possible embedding of this system in 3-dimensional Euclidean space with 
Cartesian coordinates A, Y, Z is 

1 

’71 

where x > y is 27r-periodic and sincf) = \/2x -|- 1. Here F{(l),k) is an incomplete elliptic 
integral of the 1st kind. Then = Ax{dx‘^ + dy‘^) = dX"^ + dY'^ + dZ'^. 

A general 2nd order symmetry, mod PL, can be written as Ci = 01 ^ 2 - 1 - 02^1 -|- a^lC'^, and 
the translation group generated by /C: x ^ x, y ^ y + a, leaves K? and % invariant, but 

Xi^Xi- 2an, X2^X2 + aXi - 

We classify the distinct free nondegenerate systems under this conjugacy action. We choose 
one generator Ci and determine the possibilities for C 2 such that £ 1 , C 2 , H form a quadratic 
algebra. Then we eliminate redundancies. Under conjugacy we can assume that Ci takes one 
of the 3 forms Xi + alC^, X 2 -|- a/C^, K?. 

1st case. We choose £1 = A 2 -|- aK? and try to determine the possibilities for £ 2 , up to 
conjugacy under e“^, such that £ 1 , £ 2 , form a quadratic algebra. (As we go through the cases 
step by step, we ignore systems that have already been exhibited in earlier steps.) In general 
£2 = ciAi -|- and ci, C 3 are to be determined. We must require that 

F? = hiC\ -\- b 2 C\ -\- -|- b^C\lL 2 + b^C,\C\ -|- b%C,\L 2 }i 

-|- bjC\T-i + bgC^H -|- 6g?d^£i -I- bioH?'C 2 , 



X = 2-y/xcosy, Y = 2\/xsiny, 


z=^ 

3 


for some constants 61 ,..., 610 . There are 2 possible classes: 

1. DIA: £1 = A 2 + 6/C2, £2 = Ai + i/C^, = 2iL\ + SiCiC 2 'H + AbClU + mn^Ci. This 

class is superintegrable with 


= 0 , 

G(x,y) = ln(4x). 


2 1 5x — 2b + 2iy 

X ’ 2 x{x — b + iy)' 

D = -^(x -b + iy). 


The potential of the superintegrable system is 


V{x,y) 


bi{2x -2b + iy) 62 , ^3 , , 

xYx — b + ly xYx — b + ly x 


R22 


-2>i 

X — h + iy' 


(This is missing in the tabulation in [20], but pointed out in [19] and [24].) 
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2. DIB: Cl = A'a, £2 = /C^, = -4£| - 16£i£2^, 


= 0 , 






R 22 = 


G(x,y) = ln(4x), 


D = -l 
2 


(Listed as a superintegrable system in the tabulation in [20].) The potential of the super¬ 
integrable system is 

T// N _ , &2 , h , ^ 

^ y) — ’—^^H-1—o + 04- 


X 


X xy^ 


2nd case. We choose £1 = Xi. Then the only possibility is 
3. DIG: £1 = Xi, £2 = /C2, 7^2 = 16C2n^ 


= 0 


r 22 = 0 , 


X X 

G(x,y) =ln(4x), ^ = ~\- 

(Listed in the tabulation in [20].) The potential of the superintegrable system is 

^ 6i(a;2 + y2) 63 , 632/ , ^ 

V[x, y) = ■—^- 4 , -^-^ hi. 

X XX 


5.0.2 Free Darboux 2 systems 

The space Darboux 2 (D 2 ) has free degenerate Hamiltonian 

k .2 


D2D: 




with a single Killing vector 1C = py and a basis, {T-L^K?, Xi, X 2 } for the 4-dimensional space of 
2nd order Killing tensors. Here, 


Xi = 2xpxPy + ^ {pI - x^pI) , T 2 = 2xypxPy + 


(y2 - x^)py -I- x2 (1 - y‘^)p. 


x 2 -|- 1 


The commutation relations are 


{/C,T’i} = 2(/C2-H), {lC,X2] = Xi, {Xi,X2} = 41CX2, 

and there is the functional relation 4?^T’2 -|- X^ — 4:IC‘^X2 — 4 'R 2 = q_ xhe degenerate potential is 

V{x,y) = 3 ^^ -4 62 - 

As shown in [19] The line element ds'^ can be realized as a two-dimensional surface embedded 
in three dimensions by 

^ yVx^ + I = 

X ’ X ’ 

( 2 x^ + 5x2 _|_ 8y2W^2 _j_ I 3 

Y + Z = - - - -arcsinhx, 

8 x 8 
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in which case, 

= dX^ + dY^ - dZ^ = (dx^ + dy 2 ). 

A general 2nd order symmetry, mod Ti, can be written as £i = 01^2 + a 2 Xi + 03 / 0 ^ and the 
translation group generated by/C: x—)'X,y—)-y + a, leaves K? and % invariant, but 

Ai ^ Ai - 2a'H + 2a/C^ A 2 ^ ^2 + aXi + a^(/C^ - U). 

We classify the distinct free nondegenerate superintegrable systems under this conjugacy action. 
We choose one generator £1 and determine the possibilities for £2 such that £ 1 , £ 2 , Id form 
a quadratic algebra. Under conjugacy there are 3 possible choices: £1 = A 2 + a/C^, Xi, K?. 

1st case. We choose £1 = A 2 + aK? and try to determine the possibilities for £ 2 , up to 
conjugacy under such that £ 1 , £ 2 , form a quadratic algebra. (As we go through the cases 
step by step, we ignore systems that have already been exhibited in earlier steps.) In general 
£2 = ciAi + C 3 K? and ci, C 3 are to be determined. We must require that 

'R?' = h\C\ + 62£2 + + hiiC\C2 + b^CiC\ + bQCiC2Td 

+ bjCiTd + bsC\Td + bgTd^ Cl + bmld^ C2, 


for some constants bi,... ,bio. There are 2 possible classes: 

1. D2C: £1 = A 2 , £2 = Xi, R? = 4£i£| + 16£f?^ - 16£i'R^ 

3x2 _ ^ 


= 0 , 


R 22 = _ 


A 22 = 


6 y 


x^ + 


x(x^ + 1) ’ 

G(x,y) = In 


^12 ^ - 2 ?/^) 


x2 + 1 


X^ 


x(x^ + 1) (x^ + y2^ ’ 

, D = -x(x2 + y 2 ). 


The potential of the superintegrable system is 


V(x,y) = 


bi + 


b2 


+ 


bs 


\/x2~+~^(x2 + 1) \ y + Y^x^ + y2 y - y^x^ + y2 


+ 64 . 


(Listed as a superintegrable system in [19].) 

2. D2B: £1 = ^ 2 , £2 = /C^, 7^2 = 16£i£i “ 16£i£2H + 16£2'R^ 

3x2 _ 


= 0 , 


G(x,y) = In 


A 22 = 


x^ + 1 


x^ 


x(x2 + 1 ) ’ 

D = —xy. 




x(x2 + 1 ) ’ 


The potential of the superintegrable system is 

^(^1 + 2/') + ^ + ^) + ^4- 


r 22 = 


(Listed as a superintegrable system in [19].) 

2nd case. We choose £1 = Xi. Then the only possibility is 
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3. D2A: £i = Ai, £2 = /C^, = IQCl - 2,2Cpi + , 

3x2-1 2 


= 0 , 


G(x,y) = In 


A 22 = 


x(x2 + 1 ) ’ 




x(x2 + 1 ) ’ 


R 22 = 0 , 


X2 + 1 


X^ 


D = —X. 


The potential of the superintegrable system is 

y{x,y) = 2^ r6i(x2 + 4y2) + -^ + hy] + 64- 
x"^ +1 \ x^ J 

(Listed as a superintegrable system in [19].) 


5.0.3 Free Darboux 3 systems 

The space Darboux 3 (D 3 ) has free degenerate Hamiltonian 

with a single Killing vector 1C = py and a basis, {"R,/C^, Ti, T 2 } for the 4-dimensional space of 
2nd order Killing tensors. Here, 

^ 1 X ■ 1 2 le^(e^ + 2 ) 2 

Ti = -e smy p^py + - cos y P. - ^ ^ ^ cos y Py, 

^ 1 ^ 1 e 2 ^ . 2 le^(e^ + 2 ) . 2 

T 2 = - -e cos y p^py + - sm y P. - ^ ^ ^ sm y Py. 


The commutation relations are 

{/c,Ti} = -T2, {/c,T2} = Ti, {Xi,X2} = ^}cn, 

and there is the functional relation X X 2 — — \K?% = 0. The degenerate potential is 


V{x,y) 


hi 
+ 1 


+ 62 - 


As shown in [19], we can embed D2> as a surface in 3D Minkowski space with coordinates X, 
K, Z in such a way as to preserve rotational symmetry. Let 


X = 2y/2e 2 -^1 -|- e ^ cos ^, 


Z = — In 
12 


Y = 2y/2e 2 -^ 1 - 1 - e~^ sin ^, 


/#(6 -h 5e"^)V3 + 2e2^ -h 5e^ -M\ ^ ^ , --- 

- _ e-^v/2^3 + 2e2^ + 5e*. 

V ^(6 + 5e*)V3 -h 2e2^ -h 5e^ - 1 / 


Then 


dX"^ + dY^ - dZ^ = ^ (dx 2 + dy 2 ) ^ 


An alternate basis is {'H,X^,yi,y 2 }, where 




(X ^ 2 ie*(e*-|-2) 2 

U = ( e pd>, + ^—^P, - 2 ^.,,1 I e 


ly 


2 + 1 
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3^2 = 


r * 2 , ie^{e^ + 2) 2\ -™ 

Pd>,-2—^P, + ^ e» + l 


and Xi = —\{yi — 3 ^ 2 )) d :2 = ~ j(3^i + 3 ^ 2 )- The new commutation relations are 
{K,,yi} = iyi, {/C,3^2} = -*3^2, {3^i, 3^2} =-4i/c?^, 

and the functional relation is 3 ^i 3^2 — — 2}C^'H = 0. 

Returning to the first basis, we note that a general 2nd order symmetry, mod Ti, can be 
written as £1 = 01^2 + fl2'3’i + as/C^. and the translation group generated by /C: x —)■ x, 
y —)• y + a, leaves and T-L invariant, but Xi —)• cosaAfi — sinaR 2 , X 2 —t sinadfi + cosaR 2 - 
We classify the distinct free superintegrable systems under this conjugacy action. We choose 
one generator Ci and determine the possibilities for £2 such that £ 1 , £ 2 , % form a quadratic 
algebra. Under conjugacy the choices are Ci = Xi + aK?, Xi + 1 X 2 -, Xi + iX 2 — K?, Xi, K?. 

1st case. We choose £1 = dfi + a/C^ and try to determine the possibilities for a, £ 2 , up to 
conjugacy under such that £ 1 , £ 2 , form a quadratic algebra. (As we go through the cases 
step by step, we ignore systems that have already been exhibited in earlier steps.) 

There are 3 possible classes: 

1st class. 

1 . D 3 A: Ci=Xi- £2 = A2 + i/C 2 . 7^2 = -£3 + i £3 _ 1^3 ^ ^£ 2 ^^ _ £^£2 + 1 (£2 + 


CDn + \Ci7e - \ll 2 }i^. Here, 

€■^( — 6* _j_ ^-2x-iy _ g-iy _|_ g-a;^ 

(1+ e-^)(e* + l)(e^+ €-* 3 ^) ' 


= 0 , 


a 22 = 


R 22 = 


3ie-^y 

ga: _g Q-iy ’ 




1 e"^(-l + 2e-^-^y + e"* + 


(l + e-*)(e^ + e-*2^) 


G(x,y) = In (2(e-" + D =-{e^ - 2e*^). 

The potential of the superintegrable system is 


U(x,y) = 


+ 


626 ^ 


+ 


b^e^+iy 


1 + e^ VrT2e^^(l + e^) \fl + 2e^ {1 + e^) 

(This is a superintegrable system listed in [19].) 

2. D3B: £1 = Xi, £2 = ^2, 7Z^ = -1?^3 + 1 (£2 ^ £2)^^ 

1 -e-"^ 


+ 64 - 


= 0 , 


A 22 = 


= 0 , 


2x 


1 + e-^ ’ 

G(x,y)=ln(2(e-" + e-2-)), D=^ 

The potential of the superintegrable system is 




2 (1 + e-^)’ 


U(x,y) = 


1 


bi + e 2 ( 62 cos I + 63 sin | 


+ &4- 


(Listed as a superintegrable system in [19].) 

3. D3C: £1 = Xi, £2 = /C2, 7Z^ = -4£2£2 + 2£l?{ + £ 2 ^^, 


A^^ = 0 , 


A 22 = 


1 — e 
1 + e“ 


^12 ^ 1 + 2 e" 


= —3coty, 


1 + e- 
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G(x,y)=ln(2(e-" + e-2-)), 

The potential of the superintegrable system is 


D = — — siny. 


V{x,y) = 


1 


6 i+e^ 


&2 h 


cos 


2 y 


■ 2 V 
sm I 


+ 64 - 


(Listed as a superintegrable system in [19].) 

2nd class. We choose Ci = Xi + 1 X 2 - Then the only new possibility is 
4. D3D: Ci = Xi + iX2, £2 = 7^2 = -ACjC2, 

1-e-^ 


= 0 , 


T 22 = 


1 + e^ 

G(x,y)=ln(2(e-" + e-2-)), 


^12 ^ 1 

1 + e- 

n — 

46 . 


5^2 = -3i, 


The potential of the superintegrable system is 

^ + h. 

(Listed as a superintegrable system in [19].) 

3rd class. We choose £1 = Ti + iT ’2 — K?. There are no new possibilities. 


5.0.4 Free Darboux 4 systems 


The spaces Darboux 4 {D 4 {b)) have free degenerate Hamiltonian 

D4(b)D: n = (5.1) 

2 cos 2 x + o' 

b 7 ^ ±2, with a single Killing vector JC = py and a basis, {"R,/C^, Ti, T 2 } for the 4-dimensional 
space of 2nd order Killing tensors. Here, 

Xi = + cos 2 xp 2 + sin 2xpxPy ), X 2 = e~‘^^[—T-L + cos 2xpy — sin 2xpxPy) ■ 

The commutation relations are 

{/C,Ti} = 2Ti, {/C,T2} =- 2 T 2 , {Ti,T2} =-S/C^-46/C?^, 


and there is the functional relation X 1 X 2 — JC^ — bJC^'H — 

h 


V{x,y) = 


2 cos 2x + b 


+ 62 - 


0. The degenerate potential is 


A general 2nd order symmetry, mod T-L, can be written as £1 = 01^2 + 02^1 + and 

the translation group generated by /C: x —)> x, y —)• y -t- a, leaves /C^ and T-L invariant, but 
Xi —>■ e~^'^Xi, X 2 —>■ e2"T’2. Also the reflection y —>■ —y, x —>■ x leaves % and K? fixed but 
Xi ^ X 2 . We classify the distinct free superintegrable systems under this conjugacy action. 
We choose one generator £1 and determine the possibilities for £2 such that £ 1 , £ 2 , Ti form 
a quadratic algebra. Under conjugacy the choices are £1 = _)_ aX 2 , Ti -|- ^2 + a/C^ (a = 0, 2, 

or a / 0,2), T 2 , /C^. 

1st case. We choose £1 = _|_ Qi;f 2 and try to determine the possibilities for £ 2 , up to 

conjugacy under such that £ 1 , £ 2 , H form a quadratic algebra. (As we go through the cases 
step by step, we ignore systems that have already been exhibited in earlier steps.) We first try 
£2 = A 2 -|- cXi. The class ac ^ 0 doesn’t yield a quadratic algebra. There are 3 other possible 
classes: 
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1. D4(b)A: £i = /C2, £2 = + ^ 1 , = -64£f + IQCiCl - GAbCjn - 

xi 2 n j ^22 2(—2 sin^ 2 x — 6 cos 2 x — 2 ) 


= 0 , 


G(x,y) = In 


sin 2 x( 2 cos 2x + 6 ) 

_ _ 6 

sin 2 x( 2 cos 2 x + 6 ) ’ 

—6 — 2 cos 2x\ „ sin 2x 


^12 4(—sin^ 2x + 6cos2x + 2) 


1 + e^y 
1 - e^y 


D = 


sin^ 2x7’ 2 

The potential of the superintegrable system is 


{e^y + e-^y). 


V{x,y) = 


sin^ 2 x 


2 cos 2 x + 6 Vsinh^ y sinh^ 2y J 2 cos 2 x + 6 
(Listed as a superintegrable system in [19].) 


62 


+ 


h 


+ 64 . 


2. D4(b)B: Cl = JC^, £2 = T 2 , TZ^ = 16£i£|, 
A 22 = 


= 0 , 


22 2 (— 2 sin^ 2 x — 6 cos 2 x — 2 ) 


^12 4(—sin^ 2x + 6cos2x + 2) 


G(x,y) = In 


sin 2 x (2 cos 2 x + 6 ) 
—6 — 2 cos 2 x 


sin 2 x( 2 cos 2x + 6 ) 

r 22 = 6 ^ 


sin 2 x 2 ^ 
2 


sin^ 2 x 

The potential of the superintegrable system is 

V(x, y) = ^ ^ + he^y\ + 64 . 

2 cos 2x + 6 \sim 2 x J 

(Listed as a superintegrable system in [19].) 

2nd case. We choose £4 = /C^ + 0 ^ 2 , £2 = Ti. Then the only possibility is a = 0, which is 
redundant. 

3rd case. We choose £1 = Ti + ^2 + a/C^, £2 = T 2 + cK?. We generate a quadratic algebra 
for the system 

3. D4(b)C: £1 = Ti + T 2 + 2K.‘^, C 2 = X 2 + /C^. V? = -166^^ + 16£?£2 - 16£i£i + 
166 £i£2'R - 166£2?^ - IQCiU'^, 

4 sin^ 2x + 26 cos 2x + 4 ^22 6 (e^^ — 1) 


A^^ = 0 , 


= - 


sin 2 x( 2 cos 2 x + 6 ) ’ 




2e^y cos 2 x + e'^y + 1 ’ 




2be^y cos^ 2x + 66 e^^ + 46 cos 2x + 46e'^^ cos 2x — 8 cos 2x + 16 cos^ x 
sin 2 x(e^ 2 ^ + 2e^y cos 2 x + l)( 2 cos 2 x + 6 ) ’ 

—cos 2 x + 16e^^ cos^ x — 4e^^ cos^ 2 x + 16e^^ cos 2 x 


+ 


G(x,y) = In 


sin 2 x(e^ 2 ^ + 2e^y cos 2 x + l)( 2 cos 2 x + 6 ) 
—6 — 2 cos 2 x 


sin^ 2 x 

The potential of the superintegrable system is 

61 


Z 1 = - (2 cos 2 x + 6 - 2 ?^ + 622 ^) . 


^( x , y ) — b+2 , b -2 

~—^— “T ^— 

sin^ X cos^ X 


Z + (1 - e‘^y)^/z Z + (1 + e 2 j/) Vz cos 


, ^3 e '^y\ , ^ 

H-- 1+^4? 


where Z = (1 — e^^)^ + 4e^^ cos^ x. (Listed as a superintegrable system in [19].) 
4th case. We choose £1 = Ti, £2 = 4 ^ 2 . We do not generate a quadratic algebra. 
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d:’i -d :’2 
2i 


5.1 Alternate free Darboux 4 systems 

There is an alternate form of D4:{b) that we shall employ. Set 

Afi + Ao 

x = iX, y = iY, j = iK.= pY, = —-—, T 2 = 

Then for b 7 ^ ±2 we have 

n, _ sinh^ 2X ^ 2 , ^2 ^ 

^“ 2 cosh 2 A + 6 ^^^+^'^^’ 

Ti = — cos(2y) (H — cosh 2A py'j — sin(2y) sinh 2A pxPY, 

T 2 = — sin(2y) (ji — cosh 2A py) + cos(2y) sinh 2A pxPY- (5.2) 

The structure equations are 

{j,yi} = - 2 ^ 2 , {^^, 3 ^ 2 } = 2 ^ 1 , {^ 1 ,^ 2 } = + 2 bjn, 

yf + y| -j^- bj^n -n^ = Q. 

We can embed 114(6) as a surface in 3D Minkowski space with Cartesian coordinates A, y, Z. 
Let 


X = y\/6 + 2 cosh 2rc, Y — Z = \/6 + 2cosh2x, 

^ y‘^y/b + 2 cosh2x f 2 cosh 4x + 2 + 46 cosh 2x + 6^ 
Y + Z = ———-h 


sinh 2x 


Vb +2 cosh 2rE(cosh Ax + 3 + 26cosh2rc) 


dx. 


Then 


ds^ = 


2 cosh 2x + 6 
sinh^ 2x 


{dx^ + dy"^) = -dX^ + dY^ - dZ^ 


The change of variable u = converts the integral into an elliptic integral in u that can be 
evaluated as a rather complicated sum of elementary functions and the elliptic integrals of types 
one, two and three. 

In terms of this alternate form the superintegrable systems can be expressed as: 


1. D4(b)A': £1 = J^, £2 = W, = 16£f - 16£i£i + 166£fl6 + 

2(—2 sinh^ 2A + 6 cosh 2A + 2) 


= 0 , 




^12 4(sinh^ 2A + 6 cosh 2A + 2) 


sinh 2 A (2 cosh 2 A + 6) 

G(x,y) = In 


r 22 = _6 


cos 2y 
sin2y 


sinh 2 A (2 cosh 2A + 6) ’ 

2. D4(b)B': £1 = J^, £2 = Ti - iy 2 , = -16£i£i, 


6 + 2 cosh 2 A 
sinh^ 2A 


= 0 , 


22 2(—2sinh^ 2A + 6cosh2A + 2) 

^ ■ 


^12 4(+sinh^ 2A + 6cosh2A + 2) 


sinh 2 A (2 cosh 2 A + 6) 


sinh 2 A (2 cosh 2A + 6) ’ 

G(x,y) = In 


= 6 L 


6 + 2 cosh 2 A 
sinh^ 2A 


3. D4(b)C': £1 = Ti - £2 = ^ 2 . 7Z^ = -^bB^ - 8£? - 8£i£i + 46(£f + CDU + 


A^^ = 0 , 


A 22 = - 


-4 sinh^ 2A + 26 cosh 2A + 4 
sinh 2A (2 cosh 2 A + 6) ’ 


R 22 = 


6 sin 2y 


— cosh 2A + cos 2y ’ 
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^12 cosh^ 2X — 66 + (86 + 16) cosh 2X cos 2y 

sinh 2X (2 cos 2y — 2 cosh 2X) (2 cosh 2X + 6 ) 
32 cosh^ 2X cos 2y + 4 cosh^ 2X — 16 cosh 2X 
sinh 2X (2 cos 2y — 2 cosh 2X ) (2 cosh 2X + 6 ) ’ 
6 + 2 cosh 2X' 


G{x,y) = In 


sinh^ 2X 


4. D4(b)D' (free degenerate): J, £2 = 3 ^ 2 , + 3^1 - - bJ^'^ - = 0. 


5.2 Summary and an example 

Erom the results of [21] and the calculations of the preceding sections we see that Theorems 1.1 
and 1.2 are valid for Darboux spaces. 

We use D2C as an example to give more details about how a nondegenerate superintegrable 
system with potential is induced from a free system: From Section 5.0.2 we have the free Darboux 
2 system 

Ci = X2, C2 = Xi, = 4:CiCl + i6Cln-i6Cin^. 


The potential equations are determined by 

3x2 _ 


= 0 , 


^22 = 


x(x2 + 1) ’ 


^12 ^ (3x4 + ^2 - 2y2) 


x(x 2 + 1 ) (x 2 + y 2 [) 
and the general solution is 


R 22 = 


6y 


V(x,y) = 


X 


61 + 


x 2 + y 


62 


2 ’ 


2 v^x 2 ~+^(a ;2 + 1 ) y y 
The induced classical system has a basis of symmetries 
..2 


+ 


+ •\/x2 + y2 -y + y^x 2 + y2 


+ 64. 


n = 


X 


(5.3) 


(5.4) 


X 2 + 


-(p+pj) + vw,), A = ^ 


x^ + 1 


(5.5) 


where 

_ 1 611/(1 - x 2 ) + h 2 {{-y + v/x 2 + y 2)2 + 1 ) - 63((j/ + V'x 2 + ^ 2)2 

2 (x2 + l) y^x 2 + y2 

^(2) _ 26iv/x2 + ^2 _ 5^((_^_|_ ^3-2 _^^2[)2 _ j) _ 63((?/+ v/x2 + j/2)2 _ l) 

4x2 + 4 

The Casimir is 


- 4£i£^ - i6Cjn + iQCin^ - 4(62 + b3)H^ + 1664 /:? + (62 + 63)£i - 

+ (86264 — 063 + 406263 + 86364)7^ + (1664 + 62)£i + 6i(—63 + 62)122 

+ f ^^ 1^2 + ^^ 1^3 + & 1&4 - 4636 I - 4626364 - 46264 ^ . (5.6) 
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We will not work out the details of the quantum case but merely note that the potential terms 
of the symmetries remain unchanged as do the 2 nd order kinetic energy terms, but there are now 
1st order terms: the Hamiltonian kinetic energy is replaced by the Laplace-Beltrami operator 
on D2 and the other generating symmetry operators are formally self-adjoint with respect to 
the D2 volume measure dx dy/{x^ + 1). 


6 Contractions of Darboux systems 

Recall that the scalar curvature of a space with metric ds^ = -|-dy^), where A(x, y) = 

^G{x,y) jg _ _g-G _j_ QyyCy Constant curvature spaces are just those for which C is 
constant; flat spaces are those for which C = 0. 

Theorem 6.1. A Darboux or Koenigs superintegrable system cannot be obtained as a geometric 
contraction of a nonzero constant curvature or flat space superintegrable system. 


Proof. Suppose we have a contraction of a nonzero constant curvature system with Hamiltonian 


n = 


pi+pi 


ph+pl 


Then there is a 1-parameter family of Hamiltonians %' = lc'(J)' ? where G'(l) = G, 


G' depends smoothly on e in the interval 0 < e < 1, and C(0) defines the metric of the target 
manifold for the contraction. Further, for e 7 ^ 0 the metric defined by G' will be a scalar multiple 
of a metric on the original constant curvature system. Thus we have C{e) = —e~^ [dx'x'G' -|- 
dyiyiG') for e 7 ^ 0, where C'(e) is nonzero and independent of x', y'. In the limit we obtain the 
constant < 11 ( 0 ), so the target manifold must be flat or of nonzero constant curvature, hence not 
a Darboux or Koenigs manifold. Similarly, if the original manifold is flat the target manifold 
must also be flat. ■ 


Now we investigate contractions involving Darboux superintegrable systems. From Theo¬ 
rem 6.1 such systems cannot be obtained as contractions of constant curvature systems. (How¬ 
ever, they are all Stackel equivalent to constant curvature space systems.) Thus we limit our¬ 
selves to the search for contractions such that the originating manifold is a Darboux space. In 
distinction to the case of constant curvature spaces as originating manifolds (where all quadratic 
algebra contractions were induced by Lie algebra contractions of e(2,C) and o(3,C)) here the 
Darboux spaces have only a one-dimensional Lie symmetry algebra, so we have to search for new 
contraction mechanisms. However, we can restrict our search to free systems and then induce 
the superintegrable system contractions automatically. 

To induce the contractions through Lie algebra methods we will employ the conformal sym¬ 
metry algebra for each free Darboux system generated by functions Q(x, y,Px,Py) = y)px + 
B{x,y)py that satisfy the relations 

{n,Q} = RQ{x,y)n, {}C,Q} = KQ{x,y)}C, JC=py. 

A straightforward calculation gives the same algebra in each case: 

gs : 'Px= Px, Vy =Py, M = XPx yPy, 

{Vx,Vy}=0, {Vx,M}=Vx, {ry,M}=ry. 

The Inonii-Wigner contractions and their geometric implementations are: 


1 ) Rxi By, eKi] 

2) ePx^ eL’yj Ad; 


X = x' + \, y = y', or X = x', y = y' + 7, 
X = ex', y = ey'. 
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Up to an isomorphism of the Lie algebra that is obtained by contraction of G3 there are no 
continuous one parametric contractions besides these two, [29]. Each of these contractions 
does induce a corresponding contraction of each free Darboux system and we have found no 
contractions of Darboux systems other than these. This approach is compatible with the use of 
generalizations of Inonii-Wigner contractions for the symmetry algebras o(3,C) and e(2,C) of 
constant curvature spaces. Indeed if we compute the conformal symmetry algebra for each free 
constant curvature system generated by functions Q{x,y,Px,Py) = A(x,y)px + B{x,y)py that 
satisfy the relations 

3 

{n, Q} = Rq{x, y)n, {Cj, Q} = Y1 y)C,, 

1=1 

where the Cj form a basis for either o(3, C) or e(2, C) we find o(3, C) in the case of the sphere 
and the affine algebra, the semidirect product of e(2, C) and the dilation AI, for flat space. One 
can show that the geometrical implementations of generalized Inonii-Wigner contractions of 
the affine algebra either agree with those of e(2, C) or the contractions cannot be implemented 
geometrically. See for example [12, Table XXIX] for the ordinary Inonii-Wigner contractions of 
the affine algebra. 


6.1 D1 contractions 

We list approaches to finding the contractions (not the contractions themselves). 

1. Let X = x' + ^, y = y' , %' = Then as e —)• 0 we have 

Jp" f 

%' = pI> + pI> , Al Px'Py' , ^2 « Py/{y'px' - x'pyi) -^ . 

Gontractions constructed from such limits would have flat space as the target manifold. 

2. Let X = x' , y = y' + R' = T-L. Then as e — >■ 0 we have 




Ax' 


Al Px'Py' - ^{pI +pI,) - 
../2 


^2 


Py'iy'Px' - X'py,) + {pI,+pI>). 


Gontractions constructed from such limits would have D1 itself as the target manifold. 
3. Let X = ex', y = ey', %' = Ae%. Then as e —^ 0 we have %' = and 


e^Xi=px'Py^ - ^{pI +pI,), 

/2 

6^2 = Py' {y'px' - x'py! ) - — {pi, +pI,), e/C = Py! . 

Gontractions constructed from such limits would have D1 again as the target manifold. 


6.2 D2 contractions 

We list approaches to finding the contractions (not the contractions themselves). 

1. Let X = x' + ^, y = y', R' = R. Then as e —>■ 0 we have 
R' = pI, + p ^,, eAi « 2px'Py ', €^^2 Ri -py. 

Gontractions constructed from such limits would have flat space as the target manifold. 
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2. Let X = ex', y = ey', %' = T-L. Then as e —> 0 we have %' = x'‘^{p^, + P^i) and 

eXi PS 2x'px'Py' + 2ypy,, eX2 = 2x'y'px'Py' + + x'‘^pI', eJC = py/. 

Contractions constructed from such limits would have the complex 2-sphere as the target 
manifold. 

3. Let X = x', y = y' + \, %' = %. Then as e —)■ 0 we have %' = +p^,), K! = /C, 

X'^X^',y ,p^i ,Py,) = Xi + ‘^^[KX -%), 

1 y':Px':Py') ~ X {kX — T-L) — Ti. 

Contractions constructed from such limits would have D2 as the target manifold. 

6.3 D3 contractions 

We list approaches to finding the contractions (not the contractions themselves). 

1. Let X = x' + In(^), y = y', %' = 2eTi. Then as e —)■ 0 we have /C = Py/ 

/ 

Ti' = e"' {pI, XpI), eXi ps ^ (sin y' p^, 

( 1 

eT 2 PS — ( - cos y p^ipy> X - sin y {p^, - p^, 

Contractions constructed from such limits would have flat space as the target manifold. 
In terms of flat space Cartesian coordinates X = r cos 9, Y = r sin 9 we have 

= ^, y' = 29, Ti'=p\+PY, eXi k, ^{p\ — py), 

eT2 PS KL'= ]^{XpY-Ypx) = ]^J■ 

2. Let X = ex', y = ey', Ti' = 4e^^. Then as e —> 0 we have Ti' = pi, + p^,/ and 

X y 

~ ^ {pI' - ^pI') , e2T2 ps - ]^Px'Py ', = Py' ■ 

Contractions constructed from such limits would have flat space as the target manifold. 

3. Let X = x', y = y' + i\ne, Ti' = Ti'. Then as e —)■ 0 we have 

y'lix',y ,Px',Py') - eXl, y 2 {x',y ,Px',Py') --y2, IC=Py>. 

Contractions constructed from such limits would have D3 as the target manifold. 

6.4 D4 contractions 

We list approaches to finding the contractions (not the contractions themselves). 

1. Let X = X' + hie,Y = Y', Ti' = ie^Ti, with b fixed. Then as e —)• 0 we have J = py and 
Ti' = e-^^'{p\,+pl,), J'=PY', 

4e^Ti ~ - cos(2y') {Ti' - 2e“^^ py') + 2 sm{2Y')e~'^^'px'Py ^ Py - pI^ 

4 e^T2 ~ — sin(2y')(^' — 2e“^^ py,) — 2cos{2Y')e~'^^ Px'Py' ~ —‘T'PxPy, 

where x, y are standard Cartesian coordinates. Contractions constructed from such limits 
would have flat space as the target manifold. 
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2. Let X = eX', V = eY', %' = with b fixed. Then as e —)■ 0 we have T-L' = 
and 

Yi-pl+n^ 2{X'^ - Y'\l, - AX'Y'px'Pv, 

€^2 ~ 2Y'py, + 2X'px'Py' ) = py’ ■ 

Contractions constructed from such limits would have the complex 2-sphere as the target 
manifold. 

3. For this case it is most convenient to take the Hamiltonian in the form (5.1). Let x = x', 
y = y' + ^\ne, Ti' = T-L, with b fixed. Then we have 

1 

^i{x',y',Px',Py') = -X[, X'^ix' ,y' ,p^i ,Pyi) = 6 ^ 2 . 

Contractions constructed from such limits would have D4{b) as the target manifold. 

The manifold 114(6) depends on a parameter, so we can extend the contractions that we 
consider by allowing the parameter to vary. We find the following new contractions: 

4. Let X = ^ -^ln(e), Y = ^, b = TL' = ^TL. Then as e —)• 0 we have TL' = 

+Py')^ = \ J = Py'-, and 

|Ti ~ -cos(y')('^' - e^'PY') - sin(T')e'^Vx'Py', 

|T2 ~ - sin(y')(^' - e^'pY') + cos{Y')e^'p\,p\,. 

Contractions constructed from such limits would have D3 as the target manifold. 


5. Let X = Y = 6 = —2 + e^, Ti' = Then as e —)• 0 we have Ti' —)• 


eY' 


^^(Px/ +Py'), J' = \J = PY', and 


X'^+l 


(y/2 _ X'^)pI, + x''\l - Y'’^)p\, 


rl4 


/ 2 / 


V2\ 


yi+n-j^--n^-2i 2 x'y'px'py' + 


+ 1 


3^2 


2y' 


‘^X'px'Py’ + j^/2 {py' ~ Px')^ ■ 


Contractions constructed from such limits would have D2 as the target manifold. 

6. Let X = eX' + Y = eY', b = — TL' = —2e‘^'H. Then as e —)• 0 we have TL' —>■ 
^{Px' +Py'), J' = tJ = PY', and 

v2 \ 


1 


1 


Y' 


yi+n - -j^ ^ ^ [pY'iY'px' - X'PY^ - —{ph +P^Y') I , 


4X' 


3^2 ~ ( PX'PY' - 


Y' 


{PX' +Py') 


Contractions constructed from such limits would have D1 as the target manifold. 

7. Using Hamiltonian (5.2) we let b = —2 + e"^, X = X', Y = Y', Ti —)• TL'. Then 
TL' = cosh^ X' {p\, +pI,), J' = J = py,, 

y'l = — cos{2Y'){T-L' — cosh2X'py/) — sin(2y') sinh2X'px'Pv'! 
y '2 = — sin(2y') (77^ — cosh2X'py/) + cos(2y') sinh 2X'px'Py'• 

Contractions constructed from such limits would have the complex 2-sphere as the target 
manifold. Expressed in terms of the symmetries of the 2-sphere we have 

M - Ji, T 2 = 2J7i 772, 77' = 773 
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6.5 Summary and examples of Darboux contractions 

From the results of [21] and the preceding sections we have an analog of Theorem 2.1 for Darboux 
spaces: 

Theorem 6.2. Every Lie algebra contraction of Q3 induces a contraction of a free geometric 
quadratic algebra Q on a Darboux space, which in turn induces a contraction of the quadratic 
algebra Q with potential. This is true for both classical and quantum Darboux algebras. 


Some examples follow. 


1. We describe how a Lie algebra contraction induces the contraction of D2C to 54, inclu¬ 
ding the potential terms. Recall for D2C we have the symmetries (5.5), potential (5.4), 
canonical equations (5.3) and Casimir (5.6). The coordinate implementation is defined by 
X = ey', y = ex'. In terms of the coordinates x', y' the canonical equations become 


y y'{x'^ + y'^) 


ai 22 


1 


o /2 /2 

,12 2.x' - y 


= - 


6x' 


(x'^ + y'^} 


The contraction is defined by 

^'i = ^1, Cf 2 = 2 '^!’ 

In the limit we find 
/' 2 


n = n', n' = -n. 


n' = y'^{ pI, + pI) +V' = si - 25i 53 + R', 
£;=5| + w(, c'^ = SiS2 + w:„ 


where 


R' = 


2x/x' 2 + y'2 


Cl + 


C 2 


+ 


C3 


+ C 4 


+ x' \/ — x' ^ 

and 61 = ^, 62 = C2, 63 = C3, 64 = C4. Here 

S\ = Px', S2 = x'px' + y'py', ^3 = ^ (R^ - y'‘^)pl, + x'y'py>, 


are the 1st order generators of the symmetry algebra for the free Hamiltonian T-L'q = 
y''^iPx' ~^Py') Poincare upper half plane: x real, y > 0. The Casimir becomes in the 

limit: 




/2 


- 4£i'£'/ + (C2 + C3)4 - l-R' + 

/ 2 

+ -—(C2 — Cs)C'2 + Tj^(c 2 + C3 -|- 4C4) = 0 . 

2 lb 


2. D4:(b)A' —)• D3E: We again give more details in our 2nd example, which involves changing 
the parameter b. For the degenerate system 


n 


sinh^( 2 x) ,9 9 ^ 

2cosh(2x) + 6(^"+^^) + 


a 

2 cosh( 2 x) -|- b' 


we can get 


n' 


1 e'^^' 

2 e^' + 1 


{pI' + pI') + 


/3 

e*' + 1 
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as a contraction case by taking 

f \ y' h ^ oj' Q oj ^ 

x = —--ln(e), ^ = 8en, “ = 

Then we have J' = 2J = p^,/, and y[ ^ —4eTi, 3^2 ~ —4eT2- 

3. D\A —)■ E9(a = —1): For system DIA we have 

Ci=X2 + bJC^, C2=Xi + b)C^, 

7^2 = 2icl + 8iCiC2n + Abcln + mn^Ci. 

We choose contraction 1) for DlA: x = x' + \, y = y', T-L = and let b = K Then as 
e —)■ 0 we find 

Cl C'l = -py'iy'px' - x'pyf), C 2 £'2 = Py'iPx' + iPy'), CL' = pI, + p^,, 

and with TZ' = {C'l, C' 2 } the Casimir becomes TZ'^ = 2iC'.^ + C'^l-L'. 

4. DIB —)■ E2: For system DIB we have 

Cl = T 2 , £2 = /C2, 7^2 = -4£3 _ lQCiC2n. 

We choose contraction 1) for DIB: x = x' + ^, y = y', CL = ■ Then as e —)• 0 we find 

Cl ^ C'l -^C' 2 , C'l = -py>{x'pyi - y'px'), C 2 = C '2 = pI^, cl' = pI, +p^/, 

and with TZ' = {£'^, C' 2 } the Casimir becomes CZ'^ = —4£2^ + AC'^CL'. 

5. Die —)■ £3': For system DlC we have Ci = Ti, £2 = CZ?■, CZ^ = V 0 C 2 CL?. We choose 
contraction 1) for DlC: x = x' + y = y'^ CL = jCL'. Then as e —)■ 0 we find 

Cl C'l = p'^,pyi , £2 = C '2 = pI, , CL' = pI, + Py,, 

and with TZ' = {C'l, C' 2 } the Casimir becomes TZ'^ = 0. 

6. DA{b)C' S7 as b ^ —2: For system D4{b)C' we have 

£1 = Ti - Cl = T 2 , 

7^2 = AbCL^ - 8£i3 - SCiCl + 46(£? + Cl)n + SCiV.'^. 

In the limit as 6 —t- —2 we have 

C[=2J^-n', C '2 = 2 JiJ 2 , 

where J7i, J 2 , Jz are the symmetries of the sphere and CL' = Ji + J 2 + ■ The limit 

system is ST. 


7 Tables/figures describing the contractions 

of the nondegenerate and degenerate Darboux systems 


We now tabulate the contractions of Darboux systems. The system on the left is the starting 
system, and an asterisk indicates that a system is degenerate 
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7.1 D1 contraction table 



1: Dl to flat space 

2: D1 to D1 

3: D1 to D1 

DIA 

CO 

me, DIA 

me, DIA 

DIB 

E2 

me 

DIB 

Die 

CO 

me 

Die 

DID* 

E5* 

DID* 

DID* 


contraction table 




1: D2 to flat space 

2: D2 to 2-sphere 

3: D2 to D2 

D2A 

CO 

51 

D2A 

D2B 

E2 

52 

D2A 

me 

CO 

54 

D2A 

D2D* 

E5* 

55* 

D2D* 

contraction table 




1: D3 to flat space 

2; D3 to flat space 

3: D3 to D3 

DSA 

E3' 

CO 

D3D 

D3B 

E3' 

CO 

D3B 

me 

El 

E2 

D3D 

DSD 

E8 

CO 

D3D 

mE* 

E3* 

Eb* 

D3E* 


7.4 D4 contraction table 



1: DA to flat space 

2: DA to 2-sphere 

3: DA to DA 

m{b)A 

El 

52 

DA{b)B 

D4(6)R 

E8 

51 

DA{b)B 

m{h)e 

E3' 

54 

DA{b)B 

D4.{b)D* 

E3* 

55* 

DA{b)D* 



A: DA to D3 

5: DA to D2 

6; DA to Dl 

7: DA to 2-sphere 

DA{b)A 

me 

D2B 

DIB 

59 

DA{b)B 

D3D 

D2A 

me 

52 

DA{b)e 

D3B 

me 

Die 

57 

DA{b)D* 

D3E* 

D2D* 

DID* 

53* 
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Figure 1. Diagram indicating the contractions of nondegenerate Darboux systems. 



Figure 2. Diagram indicating the contractions of degenerate Darboux systems. 


7.5 Free nondegenerate Darboux systems to degenerate systems 

The following are not contractions in the usual sense, but we show which nondegenerate systems 
restrict to degenerate ones. 


nondegenerate 

DIB 

me 

D2A 

D2B 

me 

DSD 

D4:{b)A 

D4:{b)B 


—)• degenerate 
DID* 
DID* 
D2D* 
D2D* 
DSE* 
D3E* 
m{b)D* 
D4:{b)D* 


7.6 Figures 

What follows are diagrams illustrating the contractions of Darboux systems. Boxes represent 
systems (blue boxes being nondegenerate and red boxes being degenerate), and arrows represent 
contractions. In these diagrams, certain contractions that do not affect the overall hierarchy 
have been omitted for aesthetic reasons, and a dotted arrow indicates that the limiting process 
inducing the contraction changes the free Hamiltonian. A diagram illustrating the nondegenerate 
Darboux systems which restrict to degenerate ones is also presented. 
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Figure 3 . Diagram indicating restrictions of nondegenerate Darboux systems to degenerate ones. 


For completeness we give the definitions of the constant curvature superintegrable systems 
that are targets of Darboux contractions. The following lists contain the defining relations for 
the free systems and the potentials of the superintegrable systems. 

Degenerate Euclidean targets: Here the coordinates are x, y and pi = Px, P2 = Pyi 
J = xp2- ypi. 

1 . E 3 : % = p\+pI, X = J, Ci= pI, £2 = P1P2, 

Casimir: —£| — £i(£i — = potential; V = a(x^ + y^). 

2 . E 5 : 'H=pI+pI, X =pi, £1 = Jpi, £2 = P2P1, 

Casimir: 5 (£2 + X^ — T-LX^) = 0, potential: V = ax. 


Nondegenerate Euclidean targets: 

1. El: £1 = J^, £2 = pI = 16£i£2('R - £ 2 ), R = a{x^ + y') + 4 + 


2. E2: 

3. £3' 

4. E8: 

5. E9: 


£1 = pI, £2 = P2J, E? = ACliT-L - £1), V = a( 4 x^ + y^) + / 3 x + 

£1 = pf, £2 = P1P2, E? = Q,V = a(x 2 + + / 3 x + yy, 

2 


£1 = £2 = (pi + ^P2)^ E^ = -16£i£|, V = 

\2 


0 


(x+iy) 


J + 7(x +y ), 

■y(x+2iy) 


A = (PI + ^P2)^ £2 = PiJ, E^ = - 2 A( 2 £r + m V = ^^+/ 3 y+ ^ 
Degenerate targets on the 2-sphere: Here A = yP3—zp2, J 2 = zpi—xps, Jj, = xp2—ypi. 

1. 53: n = J^ + Ji + Ji, X = A, £1 = (Ji + ^A)^ £2 = (Ji - iJ 2 ?, 

Casimir: — 2 i({% — X^)^ — £i£2) = 0 , potential: R = ^, 

2. 55: n = J^ + Ji + Ji, X = Ji+ iJ 2 , £1 = Ji, £2 = (Ji + iJ2)J3, 

Casimir: —^(£2 — 'T^£i) = 0 , potential: R = 

3. 56: n = Ji + Ji + Ji, X = J 3 , £1 = J 3 J 1 , £2 = J 3 J 2 , 

Casimir: — 4 £? + A + X‘^(X'^ ~E)) = 0 , potential: R = .'J , 

^ yx‘^+y‘^ 

Nondegenerate targets on the 2-sphere: 

1. 51: £1 = (Ji + iJ2)J3. £2 = (Ji + iJ 2 ?, E^ = -4£|, R = 

2. 52: £1 = (Ji + iJ 2 f, £2 = Ji, E^ = -16£?£2, R = 

3. 54: £1 = Ji, £2 = (Ji + iJ2)J3. = -4£i£2, R = 


(x+iy)'^ (T+d/R (x+iy)'^ ’ 


0z 


+ 


7(1-42^ 


+ 


(x+iy)'^ ~ (x+iy)^ ’ 


+ 


0z 


+ 


(x+iy)'^ {x-\-iy)Jx'^+y'^ ’ 


4. 57: £1 = Ji, £2 = J1J3, E^ = -4£3 - 4£2£i + ACpi, R = 


+ 


y'^Jx'^+y 


0^ _|_ J_ 


R’ 


£1 = Ji, £2 = Ji, E^ = -16£f£2 - 16£i£i + 16£i£2H, R = ^ + 4 + ^, 


5. 59: 
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8 Conclusions and discussion 

This paper is part of a series on 2D quadratic algebras, their classification, structure, repre¬ 
sentations, and especially, contractions as they relate to 2nd order 2D superintegrable systems. 
Of special interest are contractions that correspond to geometrical pointwise limiting processes 
in the physical superintegrable systems. As shown in [22], one of the consequences of contrac¬ 
ting between superintegrable systems is a series of limiting relations between special functions 
associated with the superintegrable systems, a special case of which is the Askey scheme for 
hypergeometric orthogonal polynomials. In [21] we studied quadratic algebras related to 2nd 
order superintegrable systems on constant curvature spaces and showed that there is a one-to- 
one correspondence between conjugacy classes of quadratic algebras in the enveloping algebras 
of e(2,C) and o(3,C), and isomorphism classes of superintegrable systems with potential. Fur¬ 
ther, we showed for constant curvature spaces that generalizations of Inonii-Wigner Lie algebra 
contractions of e(2, C) and o(3,C), induce quadratic algebra contractions that correspond to 
geometrical pointwise limiting processes in the physical systems. The procedure is rigid and 
deterministic. The present paper extends these results and shows that Darboux superintegrable 
systems are also characterized by free quadratic algebras contained in the symmetry algebras of 
these spaces and that their contractions are also induced by Inonii-Wigner contractions. Thus 
our basic results hold for all 2nd order 2D superintegrable systems. In follow-up papers, in 
preparation, we will classify abstract quadratic algebras and their contractions, including those 
not induced from Lie algebras, and study which of these relate to superintegrable systems. 

We intend to conclude this series by relating contractions of 2nd order superintegrable systems 
to limiting processes for R-separable coordinate systems for wave equations, introduced by 
Bocher in his famous 1894 thesis [1]. We will show that in 2D, Bocher’s limiting processes 
for cyclides and ellipses induce generalizations of Inonii-Wigner contractions of the so(4, C) 
conformal symmetry algebra of the 2D wave equation with potential and that these contractions 
explain the full contraction pattern for 2nd order superintegrable systems. Bocher’s limits, which 
we term Bocher contractions, apply to all dimensions n > 2, so this should provide a useful guide 
to the analysis of 2nd order superintegrable systems in higher dimensions. 
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